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, Abstract. In the paper we have developed a theory of stabihty preserving structural transfor- 

o . 

^SJ , mations of systems of second-order ordinary differential equations (ODEs), i.e., the transforma- 

^ , tions which preserve the property of Lyapunov stability. The main Theorem proved in the paper 

■ can be viewed as an analogous of the Erugin's theorem for the systems of second-order ODEs. 
Q>^ ■ The Theorem allowed us to generalize the 3-rd and 4-th Kelvin - Tait - Chetayev theorems. 

■ The obtained theoretical results were successfully applied to the stability investigation of the 
' rotary motion of a rigid body suspended on a string. 

: 

^ : 

, 1. Introduction 

^ ■ In the present paper we concentrate our attention on the following system of second-order 

I ordinary differential equations (ODEs): 

00 ; 

^. : (1.1) J(t)x(t) + (D(t) + G(t))x(t) + (P(t) + n(t))x(t) = F(t,x(t)) 

o ■ 

CN I where x [t) = col [xi [t) , X2 [t) , . . . , Xm [t)] is an unknown vector-function. It is well known 
. that a great number of dynamical systems can be approximately described by system (11. ip . 

^ ■ From the physical point of view the matrix J (t) = (t) > (the upper index T denotes 
. the operation of transposition) describes the inertia characteristics of a dynamical system; the 
matrices D (t) = (t) ,G{t) = -G^ (t) , P (t) = -P'^ (t) and H (t) = (t) represent a 
dissipative, gyroscopic, non-conservative positional and potential forces respectively; the vector- 
function F{t,x{t)) represents an external forces acting on the system. We assume that all 
matrices are of dimension m x m and their elements are real-valued functions of t G [0, oo) . 
Also, we assume that the vector-function F {t, x [t)) satisfies the condition \\F {t, x) \\ = 0(||a:|p). 
System (11.11) also will be referenced to as the second-order matrix differential equation. 
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As it was pointed out in [9], when the gyroscopic terms in system ( II .ip are periodic in t with 
some period r > then the formal apphcation of the averaging method to the system could 
result in the discarding of gyroscopic structures though this structures had some stabilizing 
effect on the system before averaging. Thus, it is very desirable to have a theoretical framework 
which allows us to transform the initial system into the system possessing the same stability 
properties and containing no gyroscopic structures. 

A similar problem can be stated regarding the non-conservative positional structures P (t) . 
An elimination of the non-conservative positional structures from system (11. ip without changing 
its stability properties has a particular interest in the case when system (11. ip is autonomous, 
that is, 

(1.2) Jx{t) + {D + G) x{t) + (P + n) x{t) = 0^ 

where x = col [xi (t) , X2 (t) , . . . , Xm (i)] is an unknown vector; here again the matrix J = > 
describes the inertia characteristics of the dynamical system and matrices D = , G = 
—G^, n = n^, P = —P^ represent a dissipative, gyroscopic, non-conservative positional and 
potential forces respectively. The all matrix coefficients in equation (II. 2p are assumed to be 
constant real matrices of dimension m x m. 

Suppose that the matrix L{t) is a Lyapunov matrix (see [5j, p. 117]). In several cases the 
application of a substitution 

(1.3) x{t) = J-^2L{t)i{t) 

to equation (II. 2p could lead us to the autonomous equatioEQ 

(1.4) l[^t) + V,i{t) + W,i{t) = Q 

which do not contain a gyroscopic (V = V^) and/or nonconservative positional structures 
[W = VF^). Since the matrix L{t) is a Lyapunov matrix, the null solutions of systems (II. 2p and 
(ll.4p are stable, asymptotically stable or unstable simultaneously. On the other hand, because 
of the symmetrical properties of system (II. 4p the stability investigation of its null solution can 
be an easier task than the stability investigation of the null solution of system (II. 2p . 

The approach to the stability investigation of the second-order system of ODEs (ll.2p which 
consists of reducing the initial problem to the problem of stability investigation of the corre- 
sponding equivalent (in the sense of Lyapunov, see |5l p. 118]) symmetric system (II. 4p . was 



^^For simplisity the vector function of external forces was not taken into account. 

^^In what follows, saying that equation (|1.4p is autonomous we also will bear in mind that the elements of 
the matrices Vi, Wi are real. 



STABILITY PRESERVING STRUCTURAL TRANSFORMATIONS 



3 



first suggested by D. L. Mingori in [19]. He has shown that such approach can be very use- 
ful and fruitful for the stabihty investigation in analytical mechanics. In [19] the author has 
considered the case when D > only. Though in [20] the results of D. L. Mingori were ex- 
tended on the case when D > the necessary and sufficient conditions providing that a given 
non-symmetric second-order system of ODEs is equivalent in the sense of Lyapunov to some 
symmetric second-order system remains unknown: both papers [T9] and pO] contain sufficient 
conditions only. 

Later the necessary and sufficient conditions providing that the autonomous system (11.21) 
can be reduced to some other autonomous system (11.41) with Wi = via substitution (II. 3p 
were found in papers [TT] and [12]. However, the results of that papers were obtained under 
additional assumptions that 



where H denotes a positive numerical parameter. 

In some cases the parameter H can be a part of matrix H. This will be the case when equation 
(II. 2p describes a perturbed motion of a gyroscopic systems installed on the platform which 
rotates around the vertical with the angular velocity u. Using assumptions (II. 5p and assuming 
that H = H(°) + HU^"\ where matrices HW,H(-^) are independent on H, the necessary and 
sufficient conditions providing the reducibility of system (II. 2p to some other system (II. 4p with 
Wi = W-[ where obtained in [T3] . 

In the present paper without any additional assumptions we have obtained the necessary 
and sufficient conditions (in terms of the matrix coefficients) providing that a given system of 
second-order ODEs is equivalent in the sense of Lyapunov to some other system of second- 
order ODEs with symmetric matrix coefficients. We have considered both the autonomous and 
non-autonomous cases. In the case when the initial system is autonomous we require that the 
reduced system be autonomous too. 

The paper is organized as follows. 

In Section |2] we introduce the notion of the structural transformation of a system of second- 
order ODEs and give the definition of the L^-equivalent systems of second-order ODEs. Us- 
ing the notion of the Lfc-equivalence we formulate two symmetrization problems for the non- 
autonomous system of second-order ODEs: the problem of Elimination of Gyroscopic Struc- 
tures (EGS problem) and the problem of Elimination of Non-conservative Positional Structures 



G = HG 



(1.5) 



dL {t) 
dH 



0,Vt > 



D > 0,det(G) 7^ 0, 
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(ENPS problem). In the section the necessary and sufficient conditions providing the solvabihty 
of the both problems where obtained. 

In Section [3] we reformulate the EGS and ENPS problems for the case of the autonomous 
systems and introduce the notion of the L-equivalence of two autonomous systems of second- 
order ODEs. Theorem [6] proved in the section can be considered as an analogous of the Erugin's 
theorem (see [5, p. 121]) for the autonomous systems of second-order ODEs. Some useful 
consequences from Theorem [6] are stated in Section [5l Among them there is a theorem which 
generalizes the theorems of Mingori (see [19]) and Miiller (see |20]). 

In Section |5] we discuss the question of the interconnection between the notions of the L^- 
equivalence and equivalence in the sense of Lyapunov. 

In Section |6] we demonstrate how the using of structural transformations can facilitate the 
stability investigation of the null solution of the autonomous second-order system of ODEs 
describing the rotary motion of a rigid body suspended on a string. 

Section [7] contains several conclusions about the theoretical results presented in the paper. 

2. Structural transformations of the non-autonomous systems of 
second-order ordinary differential equations 

Let us consider the following system of second-order ordinary differential equations: 



that A{t) ,B (t) are square matrices of order m whose elements are continuous on [tg, +c>o) 
functions, i.e., A{t) , B (t) G Mm (C [to, +oo)) . Also we will use the notation Mm (C* [to, +oo)) , 
i = 1,2 to denote the linear spaces of square matrices of order m whose elements belong to 
the functional space C" [to, +oo) , i = 1,2, and the notation „(M) will be used to denote the 
space of constant real matrices of dimension m x n. 

Definition 1. The structural transformation of the second-order system of ordinary differential 
equations (12. ip is the transformation of unknown vector x which can be expressed in the form 



(2.1) 



x + A(t)x + 5(t)x = 0, 



where x = "S" (t) = [xi {t) , . . . 



Xm (t)] is an unknown vector-function. By default, we assume 



(2.2) 



x = L(t)e, 



where ^ = [^i (t) , . . . , ^m (i)]^ is a new unknown vector-function, L (t) G (C^ [to, +oo)) , 
det (L(t)) 7^ 0, Vt G [to,+oo) . 
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Applying transformation (12.21) to system (12.11) we obtain the following system of second-order 
ordinary differential equations: 

L{t)at)+{2L{t)+A{t)L{t)^i{t) + 

(2.3) 

+ (t) + A{t)L (t) + B{t)L (t)) e (t) = 0, 
or , in more convenient form, 

(2.4) ■^+v{t)i + W{t)^ = 0, 

where 

V {t) = L-i (t) (2L (t) + A{t)L (t)) , 

(2.5) 

W (t) = L-i (t) (^L (t) +A{t)L (t) + B{t)L (t)) . 

Apparently we have that V {t) ,W (t) G Mm (C [to, +oo)) . Therefore, applying transformation 
(12. 2p to the system with continuous on [to,+oo) matrix coefficients (12. ip . we arrive at system 
(12. 4p whose matrix coefficients are continuous on [to, +oo) too. 

Definition 2. We say that the system of second-order ODEs (12. ip is Lk-equivalent to system 
(12. 4p {k G {0,1, 2}) if there exists a matrix L (t) G (C^ [to, +c>o)) satisfying conditions 

1) |det(L(t))| > r/ > 0, Vt G [to,+oo] , 



2) sup 

t6[to,+oo) 



dp ^ ' 



< +00, Vi G 0, k, 



together with equalities (12. 5p . A matrix L{t) G (C^ [to,+oo)) which satisfy conditions 1) 



2) /or some k G {0, 1, 2} zs called an Lk-matrix. 

According to the definition given in p, p. 353], a matrix L{t) G Mm(C^[to, +oo)) which 



satisfy conditions 1), 2) for A; = 0, is called a regular on [to, +oo) matrix. Transformation (12. 2p . 
where L (t) is an L2-matrix can also be referenced to as a Lyapunov transformation of system 
of second-order ODEs (compare with the definition of a Lyapunov transformation form [5l p. 
116]). 

Let us consider the following symmetrization problems for the given system of second-order 
ODEs ([21]): 

1) the problem of Elimination of Gyroscopic Structures (EGS problem) which consists in 
finding an L^-matrix L (t) {k = 0, 1, 2) together with matrices V (t) ,W (t) G (C [to, +oo)) , 
V (t) = V'^ {t) , such that equalities (123!) hold true \/t G [to, +cx)) ; 
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2) the problem of Elimination of Non-conservative Positional Structures (ENPS problem) 
which consists in finding an L^-matrix L (t) (A; = 0, 1, 2) together with matrices V (t) ,W (t) G 
Mm (C [to, +oo)) ,W{t) = (t) , such that equahties ([23]) hold true Vt e [to, +oo) . 



If the matrices L{t) ,V (t) , W (t) mentioned in items 1) and/or 2) exist then we say that the 
EGS and/or ENPS problems for system (12. ip can be solved by means of Lfc-transformation. 

Both symmetrization problems can be stated in terms of the L^-equivalence in the following 
way: 

1) to find a system (12.41) which is Lfc-equivalent to the given system (12. ip and such that 
V{t) = V^{t) (EGS problem); 

2) to find a system (12 .40 which is Lfc-equivalent to the given system (12. ip and such that 
W{t) = W^{t) (ENPS problem). 

Let us find the necessary and sufficient conditions (in terms of matrices A{t) ,B (t)) providing 
the solvability of the EGS and/or ENPS problems for the given system (12. ip . or, in other words, 
the necessary and sufficient conditions providing that system (12.10 is L^-equivalent to some 
system ([23!) with V (t) = (t) and/or W (t) = (t) for some A; = 0, 1, 2. 

Supposing that the matrix coefficient in front of the vector-function ^ in system (12.40 is 
symmetric (i.e., there is no gyroscopic structures), we arrive at the following matrix differential 
equation with respect to the unknown Lfc-matrix L (t): 

2 (l (t) (t) - L (t) (t)) +A{t)L (t) (t) - 

(2.6) 

~L (t) (t) (t) = 0. 

Similarly to that, assuming that the matrix coefficient in front of the vector-function ^ in 
system (12.40 is symmetric (i.e., there is no non- conservative positional structures) we arrive at 
the equation 

L (t) (t) - L (t) (t) + A{t)L (t) (t) - 

(2.7) 

-L (t) (t) (t) + B{t)L (t) (t) - L (t) (t) 5^ (t) = 0. 
It is easy to verify that there exists a unique pair of matrices K{t), S{t), such that 

(2.8) LL^ = L (t) (t) = K{t) + S{t), K (t) = -K^ (t) , S (t) = 5^ (t) . 

If matrix L{t) is an L^-matrix {k = 0,1,2) then matrices K [t) and S {t) (12.80 belongs to 
Mm{C^[to, +c>o)). It is easy to see that 

(2.9) I (L (t) (t)) = L (t) it) + L (t) (t) = 2S (t) , 
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and 
(2.10) 



t 

L (t) L^{t) = 2 j S {v) du + So, L (to) (to) = = > 0. 



to 



Taking into account equalities (12.81) . f l2.10p . we can rewrite equations (12. 6p and (12.71) in the 
form of 

(2.11) 4:K{t)+A{t) (^2}siiy)diy + So^ - (2 J S (u) du + Sq^ (t) = 

and 

2k {t) + A (t) (S (t) + K (t)) - (S (t) - K (t)) (t) + 



(2.12) 



-B (t) ( / 25 (z/) du + So] - ( J 2S {u) du + So] {t) = 

.to / \to J 



respectively. The question arises: what necessary and sufficient requirements have to be imposed 
on the matrices K (t) and 5" (t) to provide the existence of an L^-matrix L (t) which satisfy 
equality (12.81) ? The answer to this question is given by the following theorem. 

Theorem 1. A regular on [to, +00) matrix L (t) which satisfies equality (12.81) exists if and only 
if the matrices K {t) , S (t) belong to {C [to, +00)) and satisfy the following inequalities: 

t 



(2.13) 



2 J Tr {S (u)) du + Tr (So) 

to 



< /i^v^ G [to, +00) 



(2.14) 



t 

det (2^5 (u) du + So^ > rj^ ,^t E [to, +00) 

to 



for some constants /i > 0, 77 > and real valued positive definite symmetric matrix So G 



Mm (M) . 



Proof. Necessity. Suppose that there exists a matrix L (t) which belongs to (C^ [to, +00)) 
and satisfies equality (12. 8p together with inequalities 



(2.15) 



\L{t)[\p<^, Vt G [to,+oo)l!)|, 



(2.16) 



|det(L(t))| > ?7, VtG[to,+oo), 



■^^Here \\A\\p denotes the Frobenius norm of matrix A, that is, \\A\\p = ^/THAA^. 
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for some constants /i > 0, > 0. It easy to see that the matrices K (t) , S (t) appearing in fl2.8p 
belong to Mm (C [to, +oo)) , and the necessity of conditions ( I2.13p . fl2.14p immediately follows 
from (12.1 op . The necessity in the theorem is proved. 

Sufficiency . Suppose that K (t) , S (t) e M„(C[to, +oo)), K{t) = -K^it), S{t) = S^{t) and 
inequalities f l2.13p . f l2.14p hold true for some constants /i > 0, 77 > and some positive definite 
symmetric matrix 5*0. Assuming that the matrix L = L{t) satisfies equality (12.80 Vt G [to, +00) 
together with the initial condition 

(2.17) i:(to) = i:o, LoL^ = So, 

we arrive at the conclusion that equality (I2.10p together with inequality (12.130 imply inequality 
(I2.15P as well as inequality (I2.14p implies inequality (12.160 . 

Let us prove that the solution L = L (t) to the Cauchy problem ( 12. Sp , ( 12.170 supplemented 
by conditions ( 12.130 . (12.140 exists and is unique on [to, T] for any arbitrary T > to. If we denote 
by Aj, 2 = 1, 2, . . . , m the ascending ordered eigenvalues of matrix 5*0, that is, < Ai < A2 < 
. . . < Am, then inequality ( 12.130 implies that 

m 

(2.18) A„ < 5^ A, = Tr (So) < /xl 

i=l 

Taking into account inequality (12.18P we can obtain from inequality (12.14p the estimate 

det (So) 
Ai = i ^> 



\ \ — ,,2(m-l) 

A2 ■ ■ ■ Am /i ^ ' 

which leads us to the inequality 

Since det(L(t)) 7^ 0, Vt G [to, +00), equality (12.80 can be rewritten in the form of 
(2.19) L = F{t,L) = {S (t) + K (t)) {L-^f . 

Now we intend to show that the matrix-valued function F (t, L) satisfies conditions of the 
Picard-Lindelof theorem (see, for example, [SI P- 8]) in the rectangle 



(2.20) ^= S it,L) eRx Mm (M) : to < t < T, 



L-Lc 



E 



< -, < 5 < 1 



Taking into account that the elements of matrix-functions 5" (t) and K (t) are continuous on 
[to, +cxd), it remains only to show that the matrix- valued function F (t, L) is Lipschitz-continuous 
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on *P f l2.20p with respect to its second argument L. This fact follows from the following inequal- 

5 

ities, which are valid for any matrices Li G Mm (K) , i = 1, 2, such that — Lo||^ < — : 
= ||Lo 1 ((Li - Lo) Lo ^ + EY'- L,' ((L2 - Lo) L,' + Ey' 

/ OD .CO 



, i=0 



j=0 



(oo ^ 



< 



E 



< 



j=l \i=l 



Ll — L2 K L2 — Li 



oils 11^0 ll^; 



l-^o^lU l|Li - L2||£;E E 11^1 ~ ^oir^;'' ||Lo - L2||-^ ^ II-^o^Ib 



i-1 



< 



i=l \j=l 



(2.21) 



< 



|Li - L2|Lk^ > 2(5*""^ = IlLi - L 

I 2iii? (1-5)' " 



2|Ie ■ 



(1 - 

In the above formula we have used the equality (see, for example, [l", p. 113]) 

00 

(A + E)-' = Y,i-'^yA\ VAgM„(M), M||<1, 
and the evident identity 

n 

= J2^'^''(^-B)B'-\ \fA,B e Mm{M.) , n = 1,2,.... 

i=l 

Using (12.211) we can estimate the norm of F [t, L) on the rectangle ^ (I2.2UI) in the following 

way: 



max 

.2 



1\T 



E \\L-Lo\\p<^ 



max IIL ^ — Lq + Lq II < 



< 



max ||L - Loll^; + ||LoML < 



[1-5) \\L-La\\E< 



k5 



(1-5)' 



+ K, 



max \\F(t,L)\\p< max \\K (t) + S (t)\\p, max (L^^)^ 

{t,L)&n ^ ie[to,T] |jL-Lo|lB<-f 



E 



q3- 



Thus, the conditions of the Picard-Lindelof theorem are satisfied and the solution of the 

r 5 

Cauchy problem (I2.17p . (I2.19P exists at least on the interval Ih = [to, h] , where h = min < T, 



K 



E 



li h = T then the theorem is proved. Otherwise, ii h < T then, applying the same reasoning 
as above to equation (12.191) with the initial condition = L {h) , we arrive at the conclusion 
that the solution to the Cauchy problem (I2.17p . (12.19^ exists at least on the interval [to,2/i] . 
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Apparently, after a finite number of iterations we will prove that the solution exists on [to,^^] • 
From the arbitrariness of T it follows that the solution to the Cauchy problem fl2.17p . fl2.19p 
exists on [to, +00) . The theorem is proved. □ 

It is not hard to verify that the matrix K {t) + S if) where K {t) = —K^ (t) , S (t) = S'^ (t) is 
bounded on [to, +00] and/or belongs to Mm{C^[tQ, +c>o)) if and only if both of the two matrices 
K{t) and S{t) are bounded on [to, +C)o) and/or belong to Mm{C'' [to, +00)). Taking this fact 
into account and using Theorem [T] we can make several conclusions stated below. 



Corollary 1. An Lk-matrix L{t) {k = 1,2) satisfying equality (12. 8 p exists if and only if 
K (t) , S (t) G Mm {C^ [to, +00)) and the following conditions hold true: 

1) there exist constants yU > 0, r/ > and matrix Sq G Mm (M) , 5*0 = S'^ > satisfying 
inequalities (^AS^i . (127141) .- 



2) sup 

ie[to,+oo) 



d' 



dP ^ ' 



+ sup 

te[io,+oo) 



d' 

—S (t) 
dp ^ ' 



< +00, Vi G 0, - 1. 
Equation (12. lip and Corollary [1] imply the following theorem. 



Theorem 2. The given system of second-order ODEs (12. ip with A{t) G Mm{C^ [to, +00)) is 
Lk-equivalent {k = 0, 1,2) to some system (12.41) with V (t) = (t) if and only if there exist 
the symmetric matrices S (t) G (C^ [to, +00)) , Sq G (K) , Sq > which define the 
skew-symmetric matrix K (t) 

t 

(2.22) AK (t) = A (t) A^ (t) - A (t) A (t) , A(t)=2 j S (ly) dv + ^o, 



to 



and satisfy conditions 

1) (I2.13p . (I2.14P for some constants /i > 0, 77 > 0; 



2) sup 

te[to,+oo) 



d' 
dt'' 



-K{t) 



+ sup 

te[to,+oD) 



It' 



-S{t) 



< +00, Vi G 0, - 1, {k^O 



From Theorem [2] we obtain the following corollary. 

Corollary 2. The given system of second- order ODEs (12. ip with A{t) G (C^ [to, +00)) is 
always Lo-equivalent to some system ((23!) with V (t) = (t) . 



From equation (I2.12p and Corollary [T] we can easily obtain the theorem which gives the 
necessary and sufficient conditions for solvability of the ENPS problem. 



■*)ln the case when fc = condition 2) should be neglected. 
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Theorem 3. The given system of second-order ODEs (12.11) is L^-equivalent {k = 0,1,2) to 
some system (12. 4p with W (t) = (t) if and only if there exist the symmetric matrices S (t) G 
Mm {C^ [to,+oo)) , 5*0 G Mm (K) , 5'o > and the skew-symmetric matrix K {t) which satisfy 
the matrix differential equation 

2k {t) + A{t)K {t) + K it) [t) + 

(2.23) +^ W S{t)-S it) it) +B{t)A (t) - A (t) (t) = 0, 

t 

A{t) = 2j S (u) du + So, 



to 



and conditions 1), 2) of Theorem^ 



It is worth to emphasize that for any initial condition K (to) = Kq = —K^ G (K) the 
solution K {t) to the matrix differential equation (I2.23P is a skew-symmetric matrix. Indeed, if 
we sum up equation (I2.23j) with the transposed equation (12.231) we obtain the Cauchy problem 

2N (t) + A{t)N {t) + N (t) A^ it) = 0, 

(2.24) 

N{t) = K it) + it) , (0) = 0. 

It is easy to see that the conditions of the Picard-Lindelof theorem for the Cauchy problem (12.241) 
are fulfilled and its solution N{t) exists and is unique on [fg, +c>o) . Therefore, the problem has 
the trivial solution only, that is, {t) = 0, \/t e [to, +oo) and K (t) = -K'^ (t) , Vt G [to, +oo) . 
Such conclusion can also be obtained from the analysis of the analytical expression for the 
general solution K (t) of equation (I2.23P (see, for example, [TTJ p. 188]). 
From Theorem [3] we can easily obtain the corollary. 



Corollary 3. The given system of second-order ODEs (12. ip is always Lo-equivalent to some 
other system ([23D with W (t) = W'^ (t) . 

Combining Theorems H] and [2] we arrive at the following one. 

Theorem 4. The given system of second-order ODEs (12. ip with A{t) G Mm{C^ [to,+oo)) is 
Lk-eqmvalent {k = 0,1,2) to some system ([23]) with V (t) = (t) , W (t) = (t) zf and 
only if there exist the symmetric matrices S (t) G M„i (C^ [to,+C)o)) , 5*0 G M^. (M) , .So > 
which define the skew- symmetric matrix K (t) (\2.22\) and satisfy conditions [l)| [2)] of Theorem 
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[2] together with equality 

A{t)M^{t) = M{t)A{t), VtG[to,+oo), 
(2.25) Id 1 t 

Condition ( 12.25^ can be obtained as a result of substitution of the matrix K (t) from equation 
f l2.23p by its expression from f l2.22p . 

Remark 1. Suppose that the conditions of at least one of the Theorems^^ or^ are fulfilled. 
Then each suitable Lk-matrix L (t) can be found as the solution to the matrix differential equation 
(12. 8p supplemented with an initial condition L(to) = Lq where Lq is an arbitrary matrix form 
Mm(M), such that LqLq = Sq. Additionally to that, the matrix coefficients of the respective 
symmetrized system (\2Ah can be found via formulas (12.50 . 

3. Structural transformations of the autonomous systems of second-order 

ordinary differential equations 

Let us consider the two systems of second-order ordinary differential equations 
(3.1) x + Ax + fix = 0, A,BeMm(M), 



(3.2) ^+vi + W^ = 0, V,WeM„ 

Definition 3. We say that the given autonomous system (13. ip is L-equivalent to system (13.20 
if there exists a regular on [0, +oo) matrix L (t) which satisfies equalitie^ 

V = (t) (2L (t) + AL (i)) , 

(3.3) 

W = L-i (t) [l (t) + AL (t) + BL (t)) , Vt G [0, +00) . 
From the first equality of (13. 3p we can easily obtain 
(3.4) L (2t) = exp {-At) C exp {Vt) , C G (M) . 

It is easy to see that if the matrix L {t) (13. 4p is regular on [0, +00) (see definition on page [5]) 
then it is an L^-matrix for A; = 0, 1, 2. Hence, we can see that the notion of the L^-equivalence 
{k = 0, 1, 2) for two autonomous systems according to definition [2] is tantamount to the notion 
of the L-equivalence according to definition [3l 



^^Without loss of generality and for the sake of simplicity, in this section we consider the segment [0, +00) 
instead of Hq, +00). 



STABILITY PRESERVING STRUCTURAL TRANSFORMATIONS 13 

In this section we consider the following symmetrization problems for the autonomous systems 
of second-order ODEs (13.11) : 

1) to find an autonomous system (13.21) with V = V'^ which is L-equivalent to the given 
system (13. ip (compare with the EGS problem); 

2) to find an autonomous system (13. 2p with W = which is L-equivalent to the given 
system (13. ip (compare with the ENPS problem). 

Let us find the necessary and sufficient requirements which have to be imposed on matrices 
A, B to provide the solvability of the EGS and/or ENPS problems for autonomous system (13. ip . 

To proceed with this task we have to introduce several convenient notations. We will use the 
notation [A, B] to describe a commutator of two square matrices A and S, that is, 

[A, B] = AB- BA. 

Also, we will use the notation {v4iy42 . . . An} to describe a superposition of commutators, that 
is, 

{A.A^} = [A, Bl {A^A^ ... A J = [{A^A2 . . . A„]. 
It is easy to ensure that the commutators obey the following properties: 

(3.5) [AB,C] = {A,C]B, VA, 5, C G M„(R) : [5, C] = 0, 

(3.6) [[A, B],C] = [A, [B, C]] , VA, B,C E M^(M) : [A, C] = 0. 
It is well known that every matrix A G Mm(M) can be expressed in the form of 

(3.7) A = Ta diag [Ai {A) E^^^^ + H'^P'\ A, (A) E^^^^ + H^p^^] 
where 

(3.8) Afc (A) = at {A) + i {A) , at {A) , {A) e M, 

= 1, 2, . . . , r. Here denotes the identity matrix; all the elements of square matrix 

are zero except those in the first superdiagonal which are equal to 1. The orders of square 
matrices E^P"^ and H^p^^ are equal to the power pk of the A;-th elementary devisor of matrix A 
and T4 denotes some nonsingular matrix from Mm(ffi) (see, for example, [11, p. 152]). 

According to formulas (13. 7p and (13.80 we define 

Ar = Ta diag [ai (A) E^p^^ + H^p^\ . . . , a, (A) E^p^^ + H^p^'^] Tx\ 

(3.9) 

Aj = Ta diag [t A {A) E^p^\ . . . , i (A) E^p^^] 



STABILITY PRESERVING STRUCTURAL TRANSFORMATIONS 



14 



then 



(3.10) 



A = Ar + Ai, ArAi = AiAr. 



Using the notion of real Jordan canonical form of a real matrix (see |7i p. 184]) it is not hard 
to prove that if A G M„ (M) then Ar, Aj G M„ (R) . 

Let us consider a Jordan matrix (see, for example, ^ p. 150]) 



(3.11) 



Jr = diag [Ji (Ai) 



Js (A.)] 



where Jj (Aj) denotes a Jordan block of size rrii corresponding to the eigenvalue Aj G M, i 
1, . . . , s. For definiteness we will use the assumption that 



(3.12) 



Aj > Xj, i < j, ^ 



rrii = m. 



1=1 



In the above formula mj denotes an algebraic multiplicity of the eigenvalue Aj of matrix Jr (see 
[71 p. 58]). The following lemma holds true. 

Lemma 1. Suppose that the matrix L{t) is defined by the formula 



(3.13) 



L{t) = exp (-Jflt) Qexp (JrI) , t>0, 



where Q G (M) . Matrix L{t) f l3.13p is a regular on [0, +oo) matrix if and only if the matrix 
Q possesses the following structure: 



(3.14) 



Q 



Qn Qi2 ■ ■ ■ Qis 

O21 Q22 ■ ■ ■ Q2S 



Osl Os2 ■ ■ ■ Qs 

where matrices Qij G Mrmmj {^) satisfy the conditions 



(3.15) 



det(Qjj)^0, jf\Q, 



and Oij denotes a zero-matrix of dimension rrii x ""^j, J = 1, 2, . . . , s. 

Proof. Without loss of generality, we consider the case when s = 2, that is, when the matrix Jr 
has only two different eigenvalues Ai, A2 G M, Ai > A2 of an algebraic multiplicity mi > and 
m2 > respectively, mi + m2 = m. Let us denote 



(3.16) 



Gi = Ji (0) ,G2 = J2 (0) 
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exp (jRt) = diag 



(3.17) 



mi 



r?i2 



exp {-Jut) = (exp {JrI)) ^ = diag 



i=0 i=0 
mi 

i: 

i! 



"12 



j=0 



i=0 



Necessity. Assume that the matrix L {t) (13.131) is a regular on [0, +oo) matrix. Taking into 
account formulas (I3.17P we obtain 



L (t) = exp {~jRt) 



Qn Qi2 
Q21 Q22 



exp ( J/jt) 



Lu{t) Li2{t) 

L2l{t) L22{t) 



(3.18) 



(mi i \ / m2 . \ 

(™2 , \ /mi , \ 



/m2 



L22 (t) = E E |G'^2 • 



Vj=o 



/ m2 



\i=0 



Since the matrices E ^^r^G**; j = 1,2 are nonsingular, it is easy to see that the matrix L{t) 



i=0 



(I3.13P has unbounded norm on [0,oo) unless Q21 = O21 and matrices Ln (t) , L22 (t) , whose 
elements are polynomials of t, are constant. The latter fact implies that 



(3.19) 



Ln(t) = Qii, L22{t) = Q 



22- 



Particulary, from equalities f l3.19p it follows that det (Qjj) 7^ 0, j = 1,2. Taking into account 
the equalities 



, i=0 



, 1=0 



from ( 13.18P and (I3.19P we obtain 



(3.20) 



I E ) = I E ) g,„ J = 1, 2, W > 0. 



. 1=0 



.1=0 



Equalities fl3.2Up imply that [Gj, Qjj] = 0, j = 1, 2, and we immediately arrive at the conclusion 
about necessity of commutative equalities in formula (I3.15p . The proof of the necessity is 
complete. 

Sufficiency. Assume that the matrix Q has a structure described in (13.140 . that is, 

Qii Q12 
O21 Q22 



Q 
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L{t) = exp{-Jnt) 



{3.21] 



Qii Qi2 
O21 Q22 



exp (Jnt) 



O21 Q22 



L12 (t) 



Conditions (I3.15P together with assumption fl3.12p imply that the matrix L [t) (13.211) is regular 
on [0, +00). Hence, the sufficiency is proved and the Theorem is proved. □ 

Lemma 2. Suppose that A, V, C & Mm (K) • // the matrix L{t), defined by formula 
(3.22) L{t) = exp (-At) C exp (Vt) , t > 0, 

is regular on [0, +00) then the spectra of matrices A and V has the same real part (see definition 
in [31 p. 145] ), that is, there exists a nonsingular matrix Ci G (M), such that 

Vr = C^'ArC^. 

Proof. From the commutativity of matrices Aj and Ar (13.100 it follows that 



(3.23) 



L{t) = exp {-Ajt) exp {-ARt) C exp (V^t) exp (V/t) . 



Taking into account the definitions of matrices Aj and Vi and equality ( ]3.23p we arrive at the 
conclusion that the matrix L (t) is regular on [0, +00) if and only if the matrix 

Li(t) = exp i-ARt) Cexp (F^t) 

is regular on [0, +00). On the other hand, it is easy to see that the matrix Li{t) represents the 
general solution to the matrix differential equation (supposing that C represents an arbitrary 
matrix from space M„ 



(3.24) 



^Li (f) = Li (t) Vr - ArU {t) ■ 
at 



In [H pp. 121-125] it was proved that equation (I3.24p possesses a solution Li {t) that is a 
regular on [0, +00) matrix if and only if the matrices Ar and Vr, has the same set of elementary 
devisors. It is known (see [3 p. 185]) that if the matrices Ar, Vr G (M) has the same set of 
elementary devisors then they are similar, furthermore, the similarity matrix Ci can be chosen 
from the space Mm (M) . This completes the proof of the Theorem. □ 



Lemma 3. Suppose that A, V,C,Z& Mm (M) , the matrix 
(3.25) L{t) = exp {-At) C exp (Vt) 
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is regular on [0, +00) and 

(3.26) [Z,L{t)C-^] =0,yt>0. 
Then there exists a nonsingular matrix Ci G Mm (K) , such that 

(3.27) Vr = Cr^A^Ci, 

(3.28) C-^ZC = C^^ZCi 
and the matrix 

(3.29) Li (t) = exp {-At) Ci exp (Vt) 
is a regular on [0, +00) matrix satisfying the identity 

(3.30) [Z, Li (t) Cf ^] = 0, Vt > 0. 

Proof. Suppose that the conditions of the Lemma are fulfilled. Then, according to Lemma [21 
the spectra of matrices A and V has the same real part. Thus, there exist nonsingular matrices 
Ta,Tv e Mm (M) , such that 

A = Ta{Jr + Ia)Tx\ V = Ty{jR + Iv)Ty\ 

(3.31) 

Ia = Tx'AiTa, Iv = Ty'ViTv, [Jr, Ia] = [Jr, Iv] = 0, 

where Jr is the Jordan matrix defined in fl3.1ip . 

Let us consider the matrix L{t) (13.251) . Using notation (I3.3ip . we can rewrite it in the form 

L{t) = Ta exp {-{Jr + IaY) T^'CTy exp {{Jr + Iv)t) Ty^ = 

(3.32) 

= Ta exp {-lAt) exp (- J/jt) (tx^CTv^ exp {jRt) exp {lyt) Ty'^ . 
From Lemma [1] it follows that Tj^CTy = Q, where Q G Mm (M) is the matrix defined in 

dsn. 

Formula (I3.32p leads us to the equality 

Tx'L{t)C-'TA = exp{-lAt)x 

(3.33) X diag [exp (- Ji (Ai) t) , . . . , exp (- (A,) t)] Qx 

X diag [exp ( Ji (Ai) t) , . . . , exp ( (A^) t)] exp {lyt) Q'^. 



STABILITY PRESERVING STRUCTURAL TRANSFORMATIONS 18 

From equality fl3.33p . owing to the commutation properties flS.lSp . we get 

T^^L{t)C-^TA = exp i-lAt) Qd exp (lyt) Q-^ + Ei (t) = 

(3.34) 

= Eo (t) + El (t) , 

where 

Qd = diag [Qu, ■ ■ ■,Qss] ■ 
It is easy to see that identity ( ]3.26p can be rewritten in the form of 

(3.35) [Tx'ZTa, T^'mCTA] = [T^'ZTa, Eq (t) + Ei (t)] = 0, Vt > 0. 

It is not hard to verify that the elements of matrix Eq (t) (13.341) can be expressed as linear 
combinations of functions of type 

(3.36) sin (a t) ± cos (a t) , a, e M. 

On the other hand, the elements of matrix Ei [t) fl3.34p can be expressed as linear combinations 
of functions of type 

(3.37) fPeP^ (cos {at) ± sin (at)) , p,aeR,p<0, 
p eN[j{0} , p < m. 

If the matrix Eq (t) + Ei (t) commutates with the constant matrix T^^ZTa for all t > 
(see (13.351) ) then the same remains true for each of the summands Eq (t) and Ei (t) separately. 
Indeed, assume to the contrary that there exists a value ti > 0, such that [Eq (ti) , T^'^ZTa] 7^ 0. 
It is obvious that in this case [Ei (ti) ,T^'^ZTa] 7^ 0. Taking into account the continuity of 
elements of matrices Eq (t) , Ei (t) , we obtain 

[Eo (t) , T^'ZTa] + 0, [El (t) , T-^^ZTj^ ^ 0, 

(3.38) 

Vt G [ti - (5, ti + 5] , 

for some sufficiently small positive real number b. 

It is easy to see that each element of the matrix [Eq [t] , T^ZTa\ can be expressed as a linear 
combination of functions of type fl3.36p and each element of the matrix [Ei (t) ^T^ZTa\ can 
be expressed as a linear combination of functions of type fl3.37p . Since the functions of types 
fl3.36p and fl3.37p are hnearly independent, we conclude that 

[Eo(t) + Ei(t),T^iZr4] = 



[Eo (t) , T-/ZT^ + [El (t) , T-^^ZT^ ^ 
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for some t G [ti — 6, ti + 6] . Thus, we get a contradictions to condition f l3.35p . This contradiction 
proves the incorrectness of our assumption. Therefore we proved the identity 

(3.39) [Eoit),Tx'ZTA] =0,Vt>0. 
Setting t = in (I3.39P we obtain 

(3.40) [QdQ-\T^'ZTa] =0. 

Let us construct a matrix Ci G (M) , det (Ci) 7^ satisfying equahty (13.281) . Using 
equahty (13.401) we get 

C-'ZC = C-'Ta {T^'ZTa) T^'C = Ty {Ty^C-^TA) [T^^ZTa] x 
X [Tx'CTv) Ty' = TvQ-' {Ta'ZTa) QTy' = 

(3.41) 

= TvQu' {QdQ-') {Ta'ZTa) {QQ-^') QnTy' = 

= TvQn' {Ta'ZTa) QnTy' = [TyQ-^^T^') Z {TAQoTy') . 

From equahties (13.411) it follows that the matrix Ci satisfying condition (13.281) can be chosen in 
the following way 

(3.42) Ci = TAQnTy' G M„ (M) . 

Equality (I3.27P can be obtained from the following chain of equalities 
C^'ArCi = {TAQnTy')'' {TaJrTx') TAQoTy' = 

= TvQ~p'T^' {TaJrTa^) TAQoTy' = TyJuTy' = Vr. 

Let us prove that the matrix Li (t) (I3.29P is regular on [0, +00). Taking into account equahty 
(I3.27P and executing several elementary transformations, we get 

Li (t) = exp {-At) Ci exp {Vt) = 
= exp {—Ajt) exp (— A/jt) Ci exp (V^t) exp (Vjt) = 
= exp i-Ajt) CiC^' exp {-Ant) Ci exp (V^t) exp (Vft) = 
= exp {-Alt) Ci exp {-C^^ Andt) exp (Vr^) exp (V/t) = 
= exp {-Alt) Ci exp {-Vnt) exp (Vnt) exp {Vjt) = 
= exp {-Alt) Ci exp {Vit) . 
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It is easy to see that the matrix exp (— ^4/^) Ci exp (V/t) is regular on [0, +oo). 
Now we intend to prove identity f l3.30p . Equahties 

T^^Li (t) C^^Ta = T^' (exp {-Ajt) [TAQuTy^) exp {Vit) 



X 



X {TAQoTy'^) 'ta = exp (-J^t) Qd exp (lyt) Q-^QQ^' 



= Eo (t) {QQ-r,'). 

together with commutation identities f l3.39p and (13.401) immediately lead us to the equalities 

[Z,L,it)C^'] = [T^'ZTA,Tx'L,{t)C^'TA] = 

= [Tx'ZTA,Eo{t) (Q(Qd)"')] =0, 

which are valid for all t > 0. 

The Theorem is proved. □ 

Let us denote by the set containing all the solutions of the system of linear matrix equations 

(3.43) {ZA^^^X} =0, k = 0,l,...,n, 

where Z, A E (K) are given matrices and X is the unknown square matrix of order m. 
Theorem 5. There exists a positive integer number n < m? , such that the set equalities 

(3.44) Xn = Xk, A; = n + 1, n + 2, . . . . 
hold true. 

Proof. It is not hard to verify that the set Xn can be represented in the multi-parametric matrix 
form 

" Pn 

(3.45) X^= Y^^tlfk 

where x^klj ^'^^ constant real coefficients and are the arbitrary parameters k = 1,2, . . . ,p„ 
(see, for example, [U p. 221]), < p„ < m^. 

To begin with, we prove that if for some non-negative integer n G N |J {0} the set equality 

(3.46) Xn = Xn+l, 

holds true then equalities (I3.44p hold as well. Indeed, equality (I3.46P implies that 

= [ZA^^+'^Xn+l] = [[{ZA('''^},A] ,Xn] = 

(3.47) 

= [{ZA^-)], [A,Xn]] = [ZA^^) [A,X^]} . 
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From here and below by equality of type f l3.47p we mean the equality for every element of set 
Xn. From it follows that 

(3.48) [AA'„]CA'„. 
From equality fl3.46p and inclusion (13 .48^ it follows that 

(3.49) 

= , [A,Xr]\ = {ZA^-+^) [A,X^]] = 0, 

i.e., Xn C Xn+2- On the other hand, from the definition of the set Xn it follows that Xn ^ 
'^n+i 12 Therefore, we have 

(3.50) Xn = Xn+l = Xn+2- 

Using reasoning similar to that used above and the method of mathematical induction it is not 
hard to prove that equality f l3.46p implies equalities f l3.44p . 

Now let us prove that the non-negative integer n G N |J {0} , mentioned in the Theorem, 
exists and is less then m^. For this purpose we consider the system of matrix equations 

(3.51) =0, n = 0, 1, . . . , - 1 

with respect to unknown matrix X G Mm (M). If we would show that every solution X of 
system fl3.5ip satisfies equalities 

(3.52) =0, n = m2,m^ + l,..., 

we will prove the Theorem. 

Let us consider the process of solving of system (13.510 . Suppose that n = 0. There are only 
two possible cases (see representation ( I3.45P ): 

a) 

{ZAA'o} = 0, Vcfc gR, A; = 1,...,po, 

that is, we already have found a non- negative integer n = 0, such that equalities fl3.46p hold 
true. Therefore, as it was proved above, equalities fl3.5ip . fl3.52p hold true for all X & Xq, the 
process is completed and the Theorem is proved; 
h) equality 

(3.53) [ZAXq] = 

does not hold true for all possible values of the parameters 
Cfc G M, A; = 1, 2, . . . ,po- This means that > Po > 2, because the assumption that po = 1 
or pq = immediately leads us to the equalities Xq = cqE or Z = aE, a G M respectively 
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and we arrive at the case a). Thus, equahty (I3.53P can be viewed as a system of hnear 
homogeneous equations with respect to the arbitrary parameters c^, k = 1,2, . . . ,po- Since this 
system possesses a non-zero solution, the rank ri of its matrix satisfies the two-sided inequahty 
< ri < po- If we would solve the given system we will arrive at the matrix (set) Xi ( I3.45p . In 
addition to that (see [lEl p. 40-41]) pi = pq — ri. Therefore, pi < po, that is, the number of the 
arbitrary parameters has decreased. Again, there are only two possible cases 



does not hold true for all possible values of the parameters 
Cfc G M, k = 1,2, . . . ,pi. It means that > Po > Pi + 1 > 3. Equality f l3.54p can be viewed as a 
system of linear homogeneous equations with respect to the parameters Ck, k = 1,2, . . . ,pi. 
If we would solve this new system we will arrive at the matrix (set) f l3.45p . It is obvious that 
in this case p2 < Pi, that is, the number of the arbitrary parameters has decreased again. And 
so on. 

This process could not contain more than pq < steps. The Theorem is proved. □ 
Lemma 4. Suppose that A, V, Z,C G Mm (K) and det (C) 7^ 0. Then the commutation identity 



{ZA'-^^Xi} = 0,Vcfe eR, k = l,...,pi. 



i.e., the process is completed and the Theorem is proved; 
b) equality 



(3.54) 




(3.55) 



[Z, L (t) C-^] =0, Vt > 0, 



where 



(3.56) 




holds true if and only if the infinite system of matrix equalities 



(3.57) 




holds true. 
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Proof. Necessity. The matrix-valued function L (t) (13.561) satisfies the equahties 

[A, L(") (t) C-i] = ylL(") (t) - L(") (t) C-M = 

= AL(") (t) - (t) VC-^ + L(") (t) VC-^ - L(") (t) C-^A = 

(3.58) 

= _2L("+i) (t) c-i + L(") (t) c-i (cyc-i - A) = 

= _2L("+i) (t) + 2L(") (t) C-iL(i) (0) C-i, L^") (t) = —L (t) 
V?7, G N|J {0} or, that is the same, 

(3.59) L("+^) (t) C-^ = it) C-^L^^^ (0) C-^ - ^ [A, L^") {t) C'^] . 

Suppose that the commutation identity (13.55^ holds. Let us proved that it implies the iden- 
tities 

(3.60) {ZA^'^) (L(t)C-^)} = 0, Vt > 0, Vn G NU{0}. 

In order to prove this, we will use the method of mathematical induction with respect to n. If 
n = then identity (I3.60p coincides with (13.551) . If n = 1 then form identity (I3.55p . using (I3.59P 
and the properties of commutators ( 13. 6p . we get 

= [Z, (t) C-^] = [Z, L (t) C-'L('^ (0) C-i] - 

(3.61) ^ ^ ^ 

~l [Z, [A, L it) C-']] = -\ [[Z, A] , L it) C-'] = -\{ZA {L (t) C-')} . 

Equality (I3.6ip proves identity (I3.60p with n = 1. Let us assume that identity (I3.6UI) is proved 
for n = k > 2 and let us prove it for n = k + 1. Using equality (I3.59P and the properties of 
commutators (13. 6p . from the latter assumption we obtain 

0= [{ZA('=)},L« (t)C-i] = [{ZAW},L(t)C-iL« (O)C-i]- 

(3.62) [{ZAW} , [A, L (t) C-']] = -\ [[{ZAW} , A] , L (t) C~'] = 

= -\{ZA^'^') (L(t)C-i)}. 

Therefore, according to the principle of mathematical induction, we have that identity (13.600 
holds for all n G N U {0} . 

Taking into account the arbitrariness of n G NU{0} in formula (I3.60p . we can obtain equalities 
(I3.57P via differentiation of identity (I3.60p with respect to t and subsequent substitution t = 0. 
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Sufficiency. Suppose that equalities fl3.57p hold. Let us prove that they imply identity fl3.55p . 
If n = then from fl3.57p we get 



(3.63) 



[Z, {CVC-' -A)]=2[Z, L(^) (0) C-'] = 0. 



If n = 1 then from f l3.57p . taking into account (13.60 . (I3.58P and (I3.63p . we obtain 

= [[Z, A] , {CVC-' -A)] = 2 [Z, [A, LW (o) C-']] = 

(3.64) 

= -4 [Z, L(2) (0) C-i] + 4 \Z, (0) C-^y 
Let us assume that we already have proved equahties 



-2f [Z,L^^) (0) C-i] 



(3.65) 



[Z, L(") (0) C"^] =0, n = l,2,...,A; 



for some positive integer k >2. 

From equahties (I3.57p . assumption (13. 65 p . properties of commutators (13. 6p and equality (I3.58P 

we get 

0= {ZA^^^CVC-' - A)} = [[{ZAC^-i)},^] ,{CVC-'-A)] = 
= [{zA''-'^},[A,CVC-'-A\] =2[{ZA'~'^} ,[A,L('H0)C~']] = 
= -4 [{ZA'~'^} , L(2) (0) C-i] + 4 , (L« (0) C-') 



(3.66) 



-4[[ZA(^-i)},L(2)(0)C-i] =-4[[{ZA(^-2)},A] ,L(2)(0)C-i] = 
= -4 ^ ^(0) (0) C-i]] = 8 [{Z^C^-^)} , L(3) (0) C-'] - 

-8 , L(2) (0) C-iL(i) (0) C'^] = . . . 

. = - (-2)"^ [[Z, A] , m (0) C-i] = - (-2)*^ [Z, [A, L(^) (0) C'^]] = 
= -(-2)'+' [Z,L(^+i) (O)C-i] . 



Thus, according to the principle of mathematical induction, we have that equalities (I3.65P hold 
for every non- negative integer n G N U {0} . 
From (I3.56P it follows that the matrix series 

yi^-^ (0) c~'- 



n=0 



STABILITY PRESERVING STRUCTURAL TRANSFORMATIONS 25 

is dominated by the number series 

g. (Pii.iic.c-,ir(r ,,/(i,„i,ll,,,-.|l)|y 

n=0 ■ ^ ^ 

Thus, the matrix series is uniformly convergent on [0, +00) and its sum coincides with the matrix 
L{t)C~^. This fact together with equahties f l3.65p immediately lead us to the commutation 
identity f l3.55p . This completes the proof of the Theorem. □ 

Now we are in position to prove the main theorem of the paper. It is stated below. 

Theorem 6 (An analogue of the Erugin's theorem). Suppose that A, B,V,W & (M) . The 
two systems of second- order differential equations 

(3.67) x{t) + Ax{t) + B x{t) = 0, 

(3.68) at) + vi{t) + wat) = o 

are L- equivalent if and only if there exists a nonsingular matrix C G M„i (M) satisfying condi- 
tions 

(3.69) Vr = C-^ArC, 

(3.70) m = V'^ + C-^ {4B - A^) C, 

(3.71) { (4S - A") {CVC-' - A) } = 0, n = 0, 1, . . . , - 1. 

Proof. Sufficiency. Suppose that for some nonsingular matrix C E (M) conditions fl3.69p . 
(13.701) and f l3.7ip are fulfilled. It is easy to see that the matrix 

(3.72) L (2t) = exp (-At) C exp (Vt) = exp {-Ajt) C exp ( V}t) 

is regular on [0,+oo). Substituting the matrix L (t) fl3.72p into the first equality of (13. 3p we 
obtain the identity 

(t) hi (t) + AL (t)) = 

(3.73) 

= (t) {~AL (t) + L{t)V + AL (t)) = V. 
From the second equality of ( 13. 3p we get 

(t) (L (t) + AL (t) + BL {t)] = 

(3.74) 

= (t) {AB - A^) L (t) + = AW. 
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Here we have taken into account that equahties (13.711) . according to Theorem [5] and Lemma HI 
are equivalent to the commutation identity 

(3.75) [4B - A\L{t)C-^] =0, Vt > 0. 

Since the regular on [0, +00) matrix L{t) (I3.72p satisfies conditions (13.31) . systems (I3.67P and 
(13.681) are L-equivalent. The sufficiency is proved. 

Necessity. Suppose that systems (13.671) and (I3.68P are L-equivalent. Then, according to the 
definition of the L-equivalence, there exists a regular on [0, +00) matrix L (t) , such that 

(3.76) L-^ (t) (2L (t) + AL (t)) = V, 

(3.77) L"^ (t) (^L (t) + AL (t) + BL (t)) = VF", Vt G [0, +00) . 
From (13.761) we obtain that 

(3.78) L (2t) = exp {-At) C exp {Vt) , 

where C G Mm (M) , det (C) 7^ 0. Then from (I3.77p . using formula (I3.78P and setting t = 0, we 
obtain equality (13.701) and commutation identity (I3.75p . 

Since the conditions of Lemma [3] are fulfilled, we can assume that the matrix C is chosen in 
such a way that identity (13.750 . equality (I3.70p and condition (I3.69P hold and in addition to 
that matrix (I3.78P is regular on [0, +00) . From identity (13.750 . according to Lemma (jl]), we get 
equalities (13. 71 p . The necessity is proved and the proof of the Theorem is completed. □ 

Remark 2. Suppose that A,V,C G (M) and det (C) 7^ 0. If for some non-negative integer 
n G NU{0} the spectrum of the matrix Zn = {(4-B — A^) A^")} is simple, i.e., all the eigenvalues 
of matrix Zn are different, then conditions (I3.7ip are equivalent to the equalities 

(3.79) [B, CVC'^ - A] = 0, 

(3.80) [A,CVC-^]=Q. 

Proof. It is almost obvious that conditions (I3.79P and (I3.80p imply conditions (13. 71 p . 

Suppose that conditions (I3.7ip are fulfilled and for some non-negative integer n the spectrum 
of matrix Z„ is simple. Then there exists a nonsingular matrix T, such that the matrix T~^ZnT 
is diagonal with pairwise different diagonal elements. Thus, (see |H, p. 221]) we have that 

(3.81) T-^ {CVC^ - A)T = diag [ai, . . . , cx^] , G M, i G l^. 
Using Theorem [5] and equalities (13.711) we obtain 

= [[Zn, A] , CVC~^ -A] = [Zn, [A, CVC-' -A]]. 
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Applying the same reasoning as above to the latter equalities we arrive at the following repre- 
sentation, which is similar to f l3.8ip : 

T-^ [A, CVC~^ -A]T = T~^ATT-^ {CVC-^ - A) T- 
-T-^ {CVC-^ - A) TT-^AT = TAT-^ diag [ai, . . . , a.m] - 
- diag [(Ti, . . . , am] TAT~^ = diag [n, . . . ,Trn] , n eR, i el, m. 



(3.82) 



It is not hard to verify that representation (13.821) implies equality [A, CVC^^ — A] = which, on 
its own account, implies equality (I3.80p . Additionally to that equality (I3.79P obviously follows 
from (13.711) . The proof is completed. □ 

Though Theorem [6] gives us the necessary and sufficient conditions providing that systems 
(13.671) and (I3.68P are equivalent (L-equivalent, to be precise), conditions (13. 69 p . (I3.70p and (I3.7ip 
of the Theorem do not possess the property of symmetry, which is one of the main properties 
of an equivalence relation. However, this is only the matter of the wording. In that form the 
theorem about L-equivalence will be useful in the further sections of the paper. Theorem E] can 
be reformulated in the "symmetric" form presented below. 

Theorem 7 (An analogue of the Erugin's theorem in the "symmetric" form). Suppose that 
A, B,V,W e Mm (ffi) . The two systems of second-order differential equations (13.670 and (13.681) 
are L-equivalent if and only if there exists a nonsingular matrix C G Mm (R) , such that 

(3.83) CVr = AnC 

(3.84) C {4W - V^) = {AB - A^) C, 

{ {AB - A^) } {CV - AC) = {CV - AC) { {4W - V^) F ^ } , 



(3.85) 



n = 0,l,2,...,m2-l. 



Proof. To prove the Theorem it is enough to show that conditions (13.690 - (I3.7ip are equivalent 
to conditions (I3.83P - (I3.85p . It is easy to see that condition (I3.69P is equivalent to condition 
(I3.83p . as well as condition (I3.70p is equivalent to condition (I3.84p . 

Taking into account (I3.84p . from equalities (I3.85P with n = we obtain the equalities 

{4B - ^2) {CVC-^ -A) = {CVC-^ -A)C {AW - V^) C'^ = 

(3.86) 

= {CVC-^ - A) (45 - ^2) 
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which lead us to condition fl3.7ip with n = 0. Multiplying equality ( I3.86P on from the left 
and on C from the right and rearranging the summands, we get 

(3.87) [C-^ {AB - A^) C, V] = [{AW ~ V^) , V] = C'^ [(4S - A^) , A] C. 

From equalities f l3.85p with n = 1, taking into account f l3.84p and f l3.87p . we obtain the equalities 

{(45 - A^) A] {CVC-^ -A) = {CVC-^ -A)C{{AW- V^) V} C'^ = 

(3.88) 

= {CVC-^ - A) {{4B - A^) V} 

which lead us to condition fl3.7ip with n = 1. Multiplying equality f l3.88p on from the left 
and on C from the right and rearranging the summands, we get 

(C-i {{AB - A^) A}C)V-V (C-i {(4S - A^) A} C) = 

= C-^{{AB-A^)A^^^}C. 
Combining the latter equality with (I3.87P we obtain 

(3.89) { {AW - V^) V^^^ }=C-'{{AB- A') A^'^ } C. 

Therefore we have proved that the first two equalities of (13.710 (with n = 0, 1) are equivalent 
to the first two equalities of (I3.85P (with r?, = 0, 1) respectively. Besides that we have proved 
the auxiliary equalities f l3.87p and f l3.89p . Let us assume that for some positive integer fc, 
2 < k < — 1 we have proved that the first k equalities of f l3.7ip are equivalent to the first k 
equalities of ( I3.85P (with n = 0,l,...,A; — 1) respectively and the auxiliary equality (similar to 



(1339])) 

(3.90) { {AW - V^) F(^) ] = C-^{ {AB - A^) A^^'^ } C 

holds. Then from equalities f l3.85p with n = k, taking into account f l3.90p . we obtain equalities 
{{AB - A^) {CVC-^ -A) = {CVC~^ - A) C {{AW - V^) V^''^ C'^ = 

(3.91) 

= {CVC-^ - A) {{AB - A^) V^^^} 

which lead us to condition fl3.7ip with ??, = A;. In addition to that, multiplying equahty ( I3.9ip 
on from the left and on C from the right and rearranging the summands, we obtain 

(C-i { {AB - A^) yl(^) } C) V - V {C-^{AB - A^) A'^^^ } C) = 

= C-i{(4S-yl2)^(fc+i)}C'. 
Combining the latter equality with assumption fl3.90p we get 

(3.92) { {AW - V^) y e^+i) }=C-'{{AB- A') A^'^+i) } C. 
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Therefore we have proved that the first A; + 1 equahties of (I3.7ip are equivalent to the first 
A; + 1 equahties of f l3.85p (with n = 0, 1, . . . , A;) respectively. Also, we have proved the auxiliary 
equality fl3.92p . According to the principle of mathematical induction we can conclude that 
equalities f l3.7ip are equivalent to equalities (I3.85p . provided that condition (I3.84p holds. This 
completes the proof of the Theorem. □ 

4. Consequences from Theorem [6] 

Below we have stated several consequences from Theorem [6] that are related to the question of 
symmetrization of the matrix differential equation (or, in other words, the system of differential 
equations) 

(4.1) Jx+ (D + G')x+ (P + n)x = 0, 

where J, D, G, F, H e M„ (M) , J = > 0, D = D^, U = , G = -G^, P = -P^. Let us 
denote 

(4.2) A = J-^D + G)J-^^, B = J-^P + U)J^^. 

Corollary 4. Suppose that there exist a symmetric matrix V G (K) arid a nonsingular 
matrix C G (K) satisfying conditions 

(4.3) [A, CVC-^] = ACVC~^ - CVC"^A = 0, 

(4.4) [B, A - CVC-^] =B{A- CVC-^) - {A- CVC-^) B = 0, 

(4.5) CVr = ArC. 

Then the autonomous equation (14. ip is L- equivalent to the autonomous equation 

(4.6) l+vi + Wi = Q, y,l^GM^(M), 



(4.7) W = ^V^ + C-^ (b - ^A^ ] C, 



containing no gyroscopic structures {V = V 

Corollary 5. Suppose that there exist matrices V,C& (IR) , det (C) ^ satisfying condi- 
tions dOD-ira and 

(4.8) - + c-'zc - C^Z^ {C~^f = 



where 



B--A^ 
4 
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Then the autonomous equation (14. ip is L-equivalent to the autonomous equation i\A.6\i . (14. 7p . 
containing no non-conservative positional structures {W = W^). 

Corollary 6. Suppose that there exist a symmetric matrix V G {^) O'^d a nonsingular 
matrix C G Mm (R) satisfying conditions (14.31) - (14. 5p together with the equality 

(4.9) C-^ZC -C'^Z'^ {C-^f = Q. 

Then the autonomous equation (14. ID is L-equivalent to the "symmetric" autonomous equation 

m. ^ {w = w^,v = v^). 

Corollary 7. If for some non-negative integer n the spectrum of the matrix 

Zn= {{4B-A^) A^")} 
is simple then the conditions of Corollaries^^ are the necessary ones (not only sufficient!). 

Combining Theorem [6] with the theorems of Kelvin - Tait - Chetayev it is not hard to prove 
the following theorem that can be viewed as a generalization of the Mingori's |19] and Miiller's 
[201 theorems. 



Theorem 8. Suppose that the matrices V,C e Mm (M) , det (C) ^ 0, V -\-V'^ > satisfy 
conditions 

(4.10) CVr = ArC, 

(4.11) { (4S - A^) {CVC^^ - a) } = 0, n = 0, 1, . . . , - 1, 

(4.12) - {V'^f + {AB - A^) C - (45 - A'f {C-'f = 0. 



If the symmetric matrix 

(4.13) W = + - ^A' j C 

is positive definite then the null solution of system (14. ip is asymptotically stable {in the sense 
of Lyapunov ) and if matrix (14.130 is nonsingular and has at least one negative eigenvalue then 
the null solution of system (14. ip is unstable {in the sense of Lyapunov) . 

It is easy to see that if P = then conditions (I4.10p - (I4.12p can be satisfied once we take 
V = A, C = E. In this case we would have that VF = 11. This means that Theorem [H] can be 
considered as a generalization of the 3-rd and 4-th theorems of Kelvin - Tait - Chetayev (see 
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It is not hard to verify that the conditions of the Mingori's [19] and Miiller's [20] theorems 
imphes conditions f l4.10p - ( I4.12p . However the following example shows that the converse of 
above proposition is not correct in general. 

Example. Assume that 



A = diag[A,A2], B 
Ai = T 



5i ^5^(2) 



(4.14) 



Bi=T 



ai 02 

bi 62 
-62 61 



T-\ Ao = T 



T-\ Bo = T 



0.3 0.4 

U a3b2 



T- 



ai 



ai 

'3 



"11" 




"10" 


1 




1 



J = diag[l,l,l,l], T 

Then in terms of matrix coefficients of equation (14.11) we have 

D = diag [Di , D2] , G = diag [d , G2] 



ai — a2 



02 
2 



02 
2 



, D2 



as - a4 



(4.15) 



02 + fli 

= -322^(2)^ G2 = -^^(2)^ ^(2) 



04 
2 



04 
2 

a4 + 03 

-1 " 

1 



n 



Hi 



65^^'^ n2 



P = — diaa 



362^(2) 30362 5'{2) 
2 ' 2ai 



bi-h 

62 



'2 ¥ 



62 + &1 



Ho 



, ii2 



6301—6203 03^2 
01 2ai 

0362 0362+6301 
2ai 01 



Both, the Mingori's [19] and Miiller's [20] theorems demand the commutativity of the matrices 
P and D. However, it is easy to verify that for the matrices P and D (I4.15P this condition is not 
fulfilled in general. Thus, we can't use the results of the mentioned theorems for the stability 
investigation of system (14.11) . (I4.15p . On the other hand, the matrices 



(4.16) 



l^ = diag[Vi,V2], C = diag[r,r], 
14 = a,E(') - (a2 - ^) S('\ V2 = a,E^^) - (a, - ^) S^^) , 
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satisfy conditions (14 .Sp - (14. 5 p of Theorem |8l according to which the matrix W (14.131) can be 
expressed in the form of 



(4.17) 



W 



010262 — + a?&i 010462 — + aibs 
wi = 2 , W2 = 2 ■ 



The conditions of Sylvester's criterion (see [T^ c. 99]), when apphed to the matrix W (I4.17p . 
lead us to the inequalities 



(4.18) 



wi > 0, W2 > 0. 



Inequalities (14.181) together with the conditions oi > 0, 03 > (providing that the matrix V 
(I4.16P is positive definite) describe the region of the asymptotical stability of the null solution 
of system (gH), (|4A5|) . 

5. On THE INTERCONNECTION BETWEEN THE NOTIONS OF THE Lfc-EQUIVALENCE AND THE 

EQUIVALENCE IN THE SENSE OF LYAPUNOV 

It is well known that systems (12. ip and (12. 4p can be rewritten in the form of 



(5.1) 

and 

(5.2) 



d 
It 



d_ 

dt 



x 




X 




= A* (t) 




X 




X 



A* (t) 



V* (t) 



, v*(t) 



O E 

-Bit) -Ait) 

O E 

-Wit) -Vit) 



respectively. Suppose that systems (15. ip and (15. 2p are connected by the transformation 



(5.3) 



X 
X 



L(t) 



i 



, hit) 



Liiit) Luit) 

L2lit) L22it) 



,t E [to, +00) , 



Lij it) G Mm (C^ [to, +00]) . It is not hard to verify that this would be the case if and only if the 
equalities 



(5.4) 



Lu it)-Lu it)Wit)-L2i it) = 0, 
L12 it) - L12 it) V it) + Lii it) - L22 it) = 0, 

B it) Lu it) + A it) L21 it) = L22 it) W it) - L21 it) 



(5.5) 



B it) L12 it) + A it) L22 it) = -L21 it) + L22 it) V it) - L22 it) 
hold true Vt G [to, +00) . 
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In accordance with the definition of the equivalence in the sense of Lyapunov of two systems 
of first-order ODEs that was given in [5l p. 118] we can introduce the same notion for the case 
of second-order systems. 



Definition 4. We say that the systems of second-order ODEs (12. ip and (12 .4^ are equivalent 
in the sense of Lyapunov if there exists a Lyapunov matrix (see definition in [5l p. 117]) L (t) 
(I5.3p satisfying conditions (15 ■4p . (15. 5p . 



Let us assume that the matrix L (t) (15. 3 p satisfies conditions 

Li2 it) = O, Lu (t) = L{t)e Mm (C^ [to, +oo)) , 



inf |det (L {t))\ > 0, sup 

te[to,+oo] te[to,+oo) 



dt'' 



Lit) 



< +00, Vfc e 0,2. 



X 




'Lit) 









X 




.Lit) 


Lit) _ 







Then from equahties (15. 4 p we immediately obtain that L21 it) = Lu it) , L22 it) = L (t) ; trans- 
formation fl5.3D reduces to the form 



(5.6) 



and represents a Lyapunov transformation (see definition in |5l p. 117]); conditions (15. 5p reduce 
to conditions (12. 5p . i.e., 

V it) = it) (2L it) + Ait)L (t)) , 

(5.7) 

W it) = it) (l it) +Ait)L it) + Bit)L it) 



Thus, we can conclude that if systems (12. ip and (12. 4p are L2-equivalent according to Definition [2] 
then they are equivalent in the sense of Lyipunov according to Definition |H However, it is almost 
obvious that the converse of above proposition is not correct in general. It is easy to see that the 
notion of the equivalence in the sense of Lyapunov includes the notions of the L2-equivalence 
(see Definition |2]) and the L-equivalence (see Definition |3]) as partial cases. Therefore, when we 
consider the possibility of using structural transformations to aid the investigation of stability 
of the null solution of system (12. ip . we inevitably arrive at the following general problems of 
symmetrization: 

1) for the given system (I2.ip . find a Lyapunov matrix L (t) (15. 3p and matrices V it) ,W it) G 
Mm iC [to, +00)) which satisfy the symmetry conditions V it) = it) and/or Wit) = 
W'^ it) together with equalities ( l5l4l) . ( 1531) Vt G [to, +00) ; 

2) for the given autonomous system flO) . i.e, Ait) = A e M^ (M) , 5 (t) = 5 g M„ (M) , 
find a Lyapunov matrix L (t) (15. 3p and matrices V it) = V E Mm (M) , W it) = W E 
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Mm (ffi) which satisfy the symmetry conditions V = V'^ and/or W = W'^ together with 
equalities (El]), (EH]) Vt G [0, +00) . 

In the case when systems (12.11) and (12.41) are autonomous, i.e., A* (t) = A* G , 
V* (t) = V* G (IR) , the necessary and sufficient conditions providing that they are equiva- 
lent in the sense of Lyapunov where found by Erugin (see the Erugin's theorem in [3l p. 145]): 
Two systems (15. ip and (15. 2 p (A* and V* are constant matrices of the same order) are equiv- 
alent in the sense of Lyapunov if and only if the matrices A* and V* have one and the same 
real part of the spectrum or , in other words, there exists a nonsingular matrix C G {^) , 
satisfying equality 

A* = cyic-\ 



Evidently, in general case, to check whether the conditions of the Erugin's theorem are fulfilled 
could be as difficult as to solve both systems (15. ip and (15. 2p directly. However, Theorems [2] - 
m M indicate that in some cases the question about equivalence in the sese of Lyapunov of two 
systems (15. ip and (15. 2p can be answered without necessity to solve them. 

Let as suppose that V (t) ,W (t) G Mm {C^ [to, +00)) . Then substituting the expressions 
for matrices L21 (t) and L22 (t) obtained from equations (15. 4p into equation (15. 5p . we get the 
following system of second-order matrix differential equations with respect to the unknown 
matrices Ln (t) , L12 {t) : 



^ (Z (t) + Z (t) V* (t)) + (Z (t) + Z (t) V* (t)) V* (t) + Bit)Z (t) 



d 

(5.8) 

+A (t) [Z (t) + Z (t) V* (t)) =0, Z (t) = [Ln (t) , (t)] • 

Thus, we arrive at the conclusion that systems (12. ip and (12. 4 p are equivalent in the sense of 
Lyapunov if and only if system (15.80 possesses a solution Z (t) satisfying conditions 



sup 

te[to,+oo) 



d^ 

dt^ ^ ' 



<+oo, fc = 0,l, inf |det(L(t))| > 0, 

te[to,+oo) 



(5.9) 



L2X it) = Lii {t)-U2 {t)W{t) 



L22 {t) = L12 (t) - L12 (t) V {t) + Ln (t) . 

The general problems of symmetrization (GPS) stated above have not been studied in this 
paper. However, on my opinion, the problem of finding necessary and sufficient conditions for 
solvability of the GPS can be interesting from both practical and theoretical points of view. This 
problem is significantly more complicated then the problem of finding necessary and sufficient 
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conditions for solvability of the EGS and/or ENPS problems (see definitions on pp. [5] and [T3|) . 
The main reason for that is the significant complexity of conditions (15. Sp . (15. 9p for finding the 
matrices L (t) , V {t) , W {t) . On the other hand, as it was mentioned above, in some cases to 
solve the GPS for the given system the one should be able to determine the Jordan canonical 
form of the system's matrix (see the conditions of the Erugin's theorem). Evidently, in this 
case the using of structural transformations can't facilitate the stability investigation of the null 
solution of the system. 

6. Application of the structural transformations to the stability 
investigation of dynamical systems 

The stability of rotary motion of a rigid body suspended on a string. Let us consider 
the symmetrization problem for the system of second-order differential equations describing the 
perturbed motion of a heavy, symmetric rigid body suspended to the stationary point O by the 
inextensible weightless string. We assume that the string attaches to the body at the point S 
lying on the body's symmetry axis. We denote the distance between point S and the center of 
mass of the body by a, and the length of the string by b. It is known that the rotary motion of 
the body can be approximately described by the following equations (see equations (2.8) and 
(2.9) from [8j): 

JiXi + Ail + cxi — (2Ji — J-i)ujX2 + (Ai — X)ujX2 + mgaxs = 0, 
J1X2 + Xx2 + CX2 + (2 Ji — Ja)^^! — (Ai — X)ujxi + mgax^ = 0, 

(6.1) < 

mlP'x^ + mb{g — hu'^)x^ — 2mb'^uJX4 + mgaxi = 0, 
mb'^Xi + mb{g — bu'^)x4 + 2mb'^ujx^ + mgax2 = 0, 

where c = mga (e + 1) + ( J3 — Ji) w^, a = be, X = fDi, Ai = fD^. In equations (16.11) by a; > 
we denote the angular velocity of rotation of the body, by m — the mass of the body, hj g — 
the free fall acceleration, and by J* = diag [Ji, Ji, J3] — the central tensor of inertia of the body. 
The authors of [8] assume that the body is effected by the dissipative moment = —fDu, 
where D = diag [Di, Di, D3] , Z^i > 0, D3 > 0, / > 0. Additionally, we assume that 2Ji — J3 7^ 0, 
Ji > 0. 

It is easy to see that system (16. ip is equivalent to system (14. ip with 

J = diag [Ji E^^\mb^ ^ ^ ^ ^i^g 0(2)] , 

G = diag [(2Ji - J^P S^^\2mb^uj S'^^^] , 

c mga E^^^ 

mga E'^^^ mb{g - bu'^) E^^) J ' 



n = 
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F = diag [-(Ai-A)a; S^^\0^^^ , 
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" 1 " 




" -1 




1 




1 



where O'^^-' denotes the square zero matrix of order 2. Furthermore, using notation (14 ■2p we get 

A = J-^ {D + G) J-5 =Di + Gi = 



(6.2) 



diag 



A -(2Ji - Jz)u 

(2Ji - J^)u: A 





" -1 " 




,2uj 


1 





(6.3) 



B = j-5 (p + n) j-5 = Pi + Hi 



b\/Jl 



gaVm rp{2) 



Let us find out the sufficient conditions in terms of the parameters of system (16. ip which 
provide that the system is equivalent to some other system does not containing the gyroscopic 
structures and (or) non-conservative positional structures. 

The elimination of the gyroscopic structures. It is easy to verify that the spectrum of 
matrix Zq = (45 — A^) (16. 2p . (16. 3p is simple. Therefore, according to Corollary [7] the conditions 
of Corollary m are necessary and sufficient simultaneously. Let us check whether the conditions 
of Corollary m are fulfilled. The matrix CViC~^ satisfying condition (14. 3 p can be expressed in 
the form of 



(6.4) 



CVC-^ = diae 









V33 V34 






_ -Vl2 VU _ 


1 


_ -V34 V33 _ 





Taking into account representation (16.40 and the fact that matrix V is real we arrive at the 
conclusion that condition (14. 5 p can be satisfied if and only if 

A 



(6.5) 



""11 



Ji 



, V33 = 0. 



(6.6) 



CVC-^ = diag 



V12 
-V12 



Taking into account (16. 4p and (16. 5p . from condition (14. 4p we can find that 

A Jf ^ V12 + ^J^Ji^ 

-V12 - ojJ^Ji^ AJf ^ 

From formula (16. 6p it follows that the matrix is a symmetric matrix if and only if 
(6.7) J3 = 0, vi2 = 0. 

Thus, the gyroscopic structures can be excluded from system (16. ip if and only if J3 = 0. 
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We can assume that the condition J3 = is satisfied if the value of the inertia moment J3 is 
fairly small in comparison with the value of 2Ji. This can be the case when the body is heavy 
and has a shape of a cylinder with a very small transverse section. 

Following to the Sommerfeld- Greenhill concept we can set A = /iJi, Ai = /iJa, where is a 
small constant coefficient depending on the environment characteristics. Returning to the case 
of a heavy cylinder with a very small transverse section, we can assume that Ai = 0. 

The elimination of the non-conservative positional structures. As it was shown 
above, conditions (14.31) - (14.51) led us to representation (16. 6p . Let us take C = E^'^\ To satisfy 
condition (14.81) we take V12 = — and according to Corollary [5l whose conditions are fulfilled, 
obtain the matrix coefficients of equation (14. 6 p 

y = diag ^ 



1 

AJi 



(6.8) 



A^ JsojA - 2a;Ai Ji 

J3wA + 2a;AiJi A^ 

-1 



2ujXi 







w 



^(2) 



J~J\ 



/mge 



A2 



^(2) 



We see that, according to Corollary [Sj the elimination of the non-conservative positional struc- 
tures is possible without any additional restrictions on the parameters of system ( 16. ip . 

Once the non-conservative positional structures are eliminated, we can try to find out the 
region of the asymptotic stability of the null solution of system (16. ip . Since systems (16. ip and 
(14. 6p . ( 16. 8p are L-equivalent, that is, equivalent in the sense of Lyapunov, their regions of the 
asymptotic stability coincide. Let us find the region of the asymptotic stabihty of system (14. 6p . 

dnn). 

Unfortunately, the matrix V (\6.8\) is not a positive definite matrix, that is, the conditions of 
Theorem[8]are not fulfilled. However, this problem can be overcame. First of all let us emphasize 
the fact that if the parameters of system (16. ip are chosen in such a way that det (W) = then the 
null solution of system ( 16. ip is unstable. Thus, we can assume that the matrix W is nonsingular. 

It is easy to verify that if 



(6.9) 



W>0 



then function V {^) = ^"^^ + ^'^WC, where ^ = ^{t) represents an arbitrary solution of system 
(14. 6p . (16. 8p . satisfies the conditions of the Krasovsky theorem on asymptotic stability (see, for 
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example, [HI p. 42]). On the other hand if the symmetric matrix W (16.81) is nonsingular and 
has at least one negative eigenvalue then the function —V (^) satisfies the conditions of the 
Krasovsky theorem on instability (see, for example, [THl p. 51]). Thus, we can conclude that 
condition (16. 9p describes the required region of asymptotic stability. 

The conditions of Sylvester's criterion, when applied to the matrix W (16. Sp . lead us to the 
following system of inequalities: 

(6.10) P>0, PS-R^>0, 

where 



b A2 ; ' ^ • 

Returning to the case of a heavy cylinder with a very small transverse section and setting 
Ai = we see that 



Therefore, the first inequality of (I6.10p is fulfilled and the second one reduces to the form 

It is worth to emphasize that conditions (I6.10p are in good agreement with the similar conditions 
obtained in [15]. On the other hand, a sophisticated method proposed in paper |8| for the 
stability investigation of the null solution of system (16. ip results in a set of inequalities which 
do not describe the region of asymptotic stability of the system (contrary to the expectations 
of the authors of the paper). The reason for that is an essential error introduced in [8j by the 
authors. 

7. Conclusions 

In the present paper we have extended and generalized the results of a series of papers devoted 
to the question of stability investigation of the null solution of systems of second-order ODEs via 
the stability preserving structural transformations. The series was started with D.L. Mingori 
[19] and then continued by Von P. C. Miiller [20], V.N. Koshlyakov [lOj, V.N. Koshlyakov and 
V.L. Makarov [H [IH [121 [13] , V.N. Koshlyakov and V.A. Storozhenko [I5]. 

In the paper we have found the necessary and sufficient conditions providing that a given 
autonomous (non-autonomous) system of second-order ODEs is equivalent in the sense of Lya- 
punov to some autonomous (non-autonomous) system of second-order ODEs which do not 
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contain gyroscopic and/or non- conservative positional structures. Theorem [8] proved in the pa- 
per generahzes the 3-rd and 4-th Kelvin - Tait - Chetayev theorems as well as the Mingori's 
[T9] and Miiller's [20] theorems. In Section |6] it was shown that the theoretical results presented 
in the paper can be successfully applied to the stability investigation of the null solution of 
systems of second-order ODEs. 
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